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ABSTRACT 

In this paper we propose a new statistic capable of detecting non-Gaussianity in the 
CMB. The statistic is defined in Fourier space, and therefore naturally separates an- 
gular scales. It consists of taking another Fourier transform, in angle, over the Fourier 
modes within a given ring of scales. Like other Fourier space statistics, our statistic 
outdoes more conventional methods when faced with combinations of Gaussian pro- 
cesses (be they noise or signal) and a non-Gaussian signal which dominates only on 
some scales. However, unlike previous efforts along these lines, our statistic is success- 
ful in recognizing multiple non-Gaussian patterns in a single field. We discuss various 
applications, in which the Gaussian component may be noise or primordial signal, 
and the non-Gaussian component may be a cosmic string map, or some geometrical 
construction mimicking, say, small scale dust maps. 
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1 INTRODUCTION 

Current theories of structure formation may be roughly di- 
vided into two classes: active and passive perturbations. Ac- 
cording to the inflationary paradigm, quantum fluctuations 
in the ver y early universe are produced during a period of 



fects through the history of the universe. The fluctuations 
in the CMB in defect models have been found to be non- 



Gaussian (Bouchet, Bennett & Stebbins 1985), even though 



the extent and strength of this non-Gaussianity is still far 
from clear. Recently a wide class of defect models have been 



inflatioir- ( 3LuiiihaidL 1990 ) and glow Lu bccunic classica l 



shown to b e in conflict with current data (Pen, Seljak & 
Turok 1997) | Albrecht, Batt ye fc Robinson 19974 |AH>rcclit 
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flation are 



still remain ([Turok 199€| ; |Durrer fc Sakellariadou 1997] |Xf 



called passive because they are seeded at some 

initial time near the Planck time and then evolve 'deter- 
ministically', or linearly . They leave their imprint on the 
CMB at last scattering (|Pccblcs 197Cj; [Hu, Sugiyama fc Silk| 



brocht et al. 1995). In fact, some of the most interestin 
( |Battye, Robinson fc Albrecht 1997j |Avelino et al. 1997 



1997]; |Hu"", 
i997| ; [Wnl 



Sugiyama 1995a ; Hu & Sugiyama 19051)1 ; |Batty^ 



hite, Scott fc Silk 1994). Although this is not strictly 



necessary, in most cases the fluctuations in the CMB tem- 
peraturi p duo to inflationary porturbat i ono form a Gauoo i& B- 



randon 



of these (based on cosmic strings) require a non-zero cos- 
mological constant of t he sort that currently favoured by 
supernova experiments (Perlmutter et al. 1998). 

Thus it is important to be able to distinguish Gaus- 
sian from non-G aussian fluctuation s . Many different tests 
are being tried (|Kogut et al. 1996 
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Luo 1994 



Magueijo 19971; |Fcrrcira, Maguci.jo fc Silk 1097 r pchmalz- 



Thbre is another class of theories of structure forma. 



ing & Gorski 1998; Winitzki & Kosowsky 1997), adapted 



Ferreira & 



tion, topologica l defects cau s ed by phase transitions in th e 
early universe (Kibble 198C; Vachaspati & Vilenkin 1984). 
Perturbatio ns caused by defects are known as active per- 
turbations (Magueijo et al. 1996; Albrecht 1996) since they 
are continually being seeded by an evolving network of de- 
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to different experimental settings, and types of signal. Our 
statistic is designed to be used for small fields where one or 
more distinctive shapes are obscured by an extra Gaussian 
component, and so are no t visible in real space. It is our ex- 
perience in previous work (Ferreira & Magueijo 1997) that in 
such situations standard statistics, based in real space, fail 
to recognize the non-Gaussianity of the signal. The idea is to 
study the statistical properties of map derivatives which are 
only sensitive to a given scale. In this way we may separate 
out different scales, some of which may have Gaussian or 
nearly Gaussian fluctuations, some very non-Gaussian. Al- 
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though this scale filtering may be achieved u s ing the wavelet 
transform ( |Ferreira, Magueijo fc Silk 1997 ; Pando & Fang 



199r j), i| n this paper we choose to use the Fourier transform 
as our scale filter. 

One problem with this approach is that the Fourier 
transform is a global transformation, and therefore can only 
recognize non-Gaussian structures globally. It may there- 
fore offer a rather contorted description for complicated net- 
works made up of essentially simple objects. We will however 
perform a transformation over the Fourier modes: a second 
Fourier transform, in angle, for modes in each ring in Fourier 
space. We will argue that by doing so we are able to recog- 
nize mostly features of individual objects, and so bypass this 
problem. This operation returns a set of quantities which are 
blind to the random orientations of individual objects. Un- 
fortunately their random positions still affect our statistic, 
so there will be a limit (albeit less stringent) on the number 
of objects in the field before our statistic becomes confused. 

The plan of this paper is as follows. In Section II we 
look at some of the ideas involved in thinking about non- 
Gaussianity. In Section III we define our statistic and look at 
some of its properties and motivations. In Section IV we ap- 
ply our statistic to practical situations in which subtle non- 
Gaussian signals are present. We consider non-Gaussian sig- 
nals corrupted by the presence of a Gaussian process, which 
can be noise or primordial signal, and which dominates on 
all but a narrow band of scales. We consider two types of 
non-Gaussian signal. We consider CMB maps obtai ned from 
string s imulati ons, implementing the algorithms in ( Smith & 
Vilenkii 1987). We also consider geometrical constructions 



mimicking, say, small scale dust maps. We show how such 
maps look very Gaussian, but fail to confuse our statistic. 



2 NON-GAUSSIAN STATISTICS 

The fluctuations in the temperature of the radiation are 
described by =jf-(n) on a two-sphere. Generally this is ex- 
p anded in spherical harmonics, but for small fields it can be 



expand 3d in Fourier modes: 



(i) 



The assumption of homogeneity and isotropy leads to 
the relation between Fourier modes 



<a*(fei)a(fe 2 )> =tf(fci -fca)P(fca) 



(2) 



where P(k) = (ja(fc)| 2 )| fe | =fc is the power spectrum and the 
averages are ensemble averages. 



Homogeneity means that 



:(x)| 2 > = <|^( a; + f,)| 2 } 



for all vectors b which means {a* (k 1 )a(k2)) = for fei 7^ kj,- 
Isotropy in real space corresponds to isotropy in Fourier 
space so the power spectrum can only depend on the mod- 
ulus of k. 

It is often assumed that the Fourier modes form a Gaus- 
sian random field. The real and imaginary parts of the 
Fourier modes of a Gaussian random field are drawn from 
a Normal distribution with mean zero and variance P(k)/2 
which means that the phases are random and the moduli 
have a \ 2 distribution. From our experience, we know that 
in real space a Gaussian random field has little structure. We 
may of course have a field that obeys equation (M but which 



clearly has structure in real space, for example a field con- 
taining squares at random positions and orientations. This 
is clearly non-Gaussian. The important point is that the av- 
erages in equation (Q) are ensemble averages, which do not 
necessarily convert to spatial or sample averages. For a non- 
Gaussian field we expect to have to take spatial averages 
over larger regions (or bigger samples) than for a Gaussian 
random field in order to get results which approximate equa- 
tion (§). 

There are various ways we can think about non- 
Gaussianity. The formal definition is the departure of the 
distribution of the modes from a Normal distribution. This 
includes cases where the phases of the Fourier modes are 
random and there is no structure in real space. If we as- 
sume the Fourier modes are independent then homogeneity 
implies that the modes have random (uniformly distributed) 
phases. By the central limit theorem the distribution rapidly 
becomes like a Gaussian as we increase the size of the field 
and spatial averages rapidly become good estimates of en- 
semble averages, i.e. the Fourier modes come close to satis- 
fying equation (0). For these cases there will be little or no 
structure in real space. 

In order to get structure in the field, we must there- 
fore introduce correlations in the modes. This means that 
spatial averages will take longer to approach ensemble aver- 
ages than for a Gaussian random field. This is because each 
realisation of the field has fewer independent modes than a 
Gaussian random field does. Therefore, to detect this type of 
non-Gaussianity we can test how well spatial averages obey 
equation (Q) compared to a Gaussian random field, rather 
than looking at the actual distribution of Fourier modes. 
Here we choose to look at a(k) in a ring of constant fc| = k 
in Fourier space, to test the isotropy of spatial averages of 
the field. 

Another reason for looking at the field in Fourier space 
is to detect non-Gaussian fluctuations superimposed on a 
Gaussian background. This happens for instance in cosmic 
string scenarios where the effects of strin gs on photons after 
last scattering (Kaiser-Stebbins effect ([Kaiser fc Stebbins 



1984)) are superimposed on the fluctuations from before last 



scattering which can be considered to be Gaussian. The field 
may be non-Gaussian on some scales but Gaussian on others, 
making it hard to see the non-Gaussianity in real space. If 
there is more power from the non-Gaussian features than 
from the Gaussian on a particular scale this should show in 
Fourier space. 



3 ANOTHER FOURIER TRANSFORM 
3.1 Definition 

We can write the Fourier coefficients as 

a{k) = p(fc)e^ (fe) (3) 

Here we will focus on the amplitudes. As we are consid- 
ering constant k we will write 



p(k) = p k (a) 



(4) 



where a is the angle of the direction of fc in polar co- 
ordinates. 

We take another Fourier transform in the angle a 
around the ring. 
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For a completely isotropic field (|a(fe)j 2 ) in one ring in 
Fourier space is constant i.e. (pi (a)) is constant in a so 
(Fk(/3)) is a 5-function in (3. 



3.2 Motivations: Repeated shapes 

This statistic is motivated by the idea of detecting patterns 
or shapes that are repeated at different places. If one par- 
ticular feature has a temperature pattern f(x) with Fourier 
transform f(k) and the field we are looking at consists of 
the feature repeated at n different positions so the field in 
real space is 



T 

i=i 

then the Fourier transform of the whole field is 

n 

a(fe) = ^/»e ife - i,i 

z=i 

and the amplitude satisfies 



p 2 (k) = 




n + 2 cos(fc ■ (pi — bj 



(6) 



(7) 



(8) 



The amplitude consists of an envelope defined by the Fourier 
transform of the individual feature together with oscillations 
which depend on the relative positions of the repeated fea- 
tures. The idea of taking the second Fourier transform is 
to get a similar envelope with oscillations in one dimension 
around the ring of constant modulus k. 

As n gets larger p 2 (k) — » n|/(fe)| 2 since the cosine av- 
erages to zero for many positions, so we can extract the 
Fourier transform of the original feature from the Fourier 
transform of the whole field. We can write this for each ring 
as pl(a) — > n\fk(a)\ 2 so we see this applies to each ring also. 

But if the shapes are repeated with different orienta- 
tions, equation (fel) becomes 



p\k) = 



E 



f(Rik] 



I 2 + 2 E [(^e/(-Rife)3?e/(i?jfe) 



bj)) 



+ Qmf(Rik)^mf(Rjk)) cos (k ■ (6, 
+ (•Stef(R i k)^mf(R j k) 

- Qmf(R i k)^ef(R j k)) sin(fc ■ (b t ~ bj))] (9) 

where the {Ri} are the matrices of the rotations. As n gets 
larger p 2 (k) — » ^\ \f(Rik)\ 2 in a similar manner to before, 
which we can write as p 2 (a) — * ^ . |/fc(a + 9i)\ 2 where 8i is 
the angle of the £th rotation. It is clear that in each ring this 
will actually tend to a constant (in a) as n becomes large if 
the angles of the rotations are uniformly distributed. This is 
merely a consequence of the increasing isotropy of the field. 

Even with just a few repeated shapes we cannot ex- 
tract the original shape as we could in the case without ro- 
tations without knowing the number and orientations of the 
repeated shapes. But unless n is very large p\(ot) should be 
distinguishable from a Gaussian. We know that these fea- 
tures will disappear when averages are taken over a large 



enough sample. The point is to get a statistic which will 
show the lack of homogeneity and isotropy on small scales 
and for small samples. 



3.3 Averaging over small fields 

Another way of using this statistic is to look at very small 
fields (i.e. those containing only a few shapes). To do statis- 
tics on these observations we want to look at many small 
fields. If we say look at several small fields each contain- 
ing one shape but at different orientations (as we would 
expect from isotropy) then we will get the same Fourier 
transform rotated each time. Mathematically it is similar 
to the calculations leading to equation ([]) except we don't 
consider translations. The real space pattern in the ith field 



is (- 



AT 
T 



(Rix) so the Fourier transform of the field 



satisfies p t (k) = p (Rik) which we write in one ring as 

Pfc(i)(«) = pl(a + Q^) (10) 
The Fourier transform taken in angle round the ring satisfies 



F m (p) = F k (p)4 
\F k(i) (p)\ 2 = \F h (/3)\ 2 



(11) 



for each field. 

Note that if we averaged over the original single- 
Fourier transformed fields {Pfc(i)(a)} before taking the sec- 
ond Fourier transform we would get 

|^F fe(j) C3)| 2 = |W)| 2 U + 2^cos(/3(fl i -6l j )) J (12) 

which tends to a constant as we add more fields. We have 
to take the second Fourier transform of each field before 
averaging to see the non-Gaussianity. This is because the 
features in p\ (a) are repeated at different a in each field and 
their sum will tend to a constant for many fields. Taking the 
second Fourier transform takes the feature from p 2 (a) and 
transforms it to a feature at the same (5 in Fk(P) for each 
field. 

The idea of looking at small fields in this way is an 
extension of the idea in (Ferreira & Magueijo 1997) of 



analysing several small fields separately. For the e xample of a 
cosmic string o n a Gaussian background used in (Ferreira & 
Magueijo 1997) the shape spectrum shows the non-Gaussian 
feature well, but the string was always considered to have the 
same orientation. If the shape spectrum were averaged over 
different directions the result would disappear. Here we can 
look at the same fields with different orientations and the 
non-Gaussianity will still be apparent. However, if we ran- 
domly choose small fields from a larger one containing many 
shapes we would have to take translations into account as 
well and the effect will not be so clear. 



4 APPLICATIONS 
4.1 Geometric shapes 

As a simple demonstration of this statistic on an obviously 
non- Gaussian field, we will consider applying it to some reg- 
ular shapes. These are intended as mock small scale dust 
maps. We expect our geometric shapes to present the same 
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1st Fourier transform 



1st Fourier transform 



6xl0 13 - 



4xl0 13 - 



2x10 





Figure 1. First Fourier transform in one ring and modulus of 
second Fourier transform for a field containing 1 square. 



Figure 2. Averages of the first and second Fourier transform 
from 100 realisations of a field containing 1 square. 



type of non-Gaussianity, and obstacles for its detection, 
the actual dust maps. 



1st Fourier transform 



First to demonstrate the ideas in Section i.3 we will look 
at fields each containing one square at different orientations. 
We first Fourier transform the whole field and take the am- 
plitude squared (p 2 (k)). For one square the amplitude is not 
affected by translation of the square. We multipl y the field 
by a G aussian window to get rid of edge effects ( Hobson & 
Maguei foTiirj )^ 



In Figure ^ we have plotted for one field p (a) and 
|-F(/9)| at k — 60 from a 200 2 array. This corresponds to 
a physical scale of k — 21600/L where L is the size of the 
real field in degrees. The y-axis is in arbitrary units. Figure ^ 
shows the same functions averaged over 100 fields. The aver- 
aged |-FQ3)| is not in fact identical to |-F(/3)| for one field, be- 
cause of numerical effects due to the discrete Fourier trans- 
form being on a lattice. The average settles down quickly, 
however: |-F(/3)| looks the same averaged over 10 fields as it 
does averaged over 100 fields. The first Fourier transform, 
p 2 (a), tends to a constant as we average over more fields as 
we would expect, but is also affected by the lattice. 

Figures g and ^ show the equivalent graphs for a Gaus- 
sian random field. Here |-F(/3)| converges rapidly and is dis- 
tinct from the result for squares. The dips in p 2 (at) are due 
to the window function superimposed in real space to avoid 
edge effects from the first Fourier transform. The Gaussian 
is treated in this way even though we generate the original 
map in Fourier space so as to model as closely as possible the 
way in which the non-Gaussian map is treated. We would 
of course expect to tend to a delta function as the 

number of fields is increased but the window effects spoil 
this. The window used is in theory circularly symmetric but 
as it is on a lattice it necessarily has some squareness in it. 
The important point is that we get different results for the 
Gaussian and non-Gaussian fields. 

Now we test the statistic on fields with many squares. 
Figures |^ and ^| show the 1st and 2nd Fourier transforms for 



3x10" 
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2nd Fourier transform 



6x10- 



4xl0 5 
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Figure 3. First and second Fourier transforms for 1 realisation 
of a Gaussian random field. 



10 and 100 squares in one field. We see that as the number 
of squares increases the result is not so distinctive. For 100 
squares the signal is barely distinguishable from that of a 
Gaussian (see Figure ^ . In fact in real space it is obviously 
not Gaussian, suggesting this is not the best way of detecting 
non-Gaussianity. Figure |^ shows the field in real space. 

The important test however, is to see if we can de- 
tect non-Gaussianity when it is superimposed on a Gaus- 
sian background so it is not visible in real space. Figure 
^| shows a real field which is the sum of a field containing 
one square and a Gaussian field. Though in real space we 
would not guess there was non-Gaussianity present, in fact 
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1st Fourier transform 



1st Fourier transform 
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2nd Fourier transform 
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Figure 4. Average of first and second Fourier transforms over 
100 realisations of a Gaussian random field. 





2nd Fourier transform 
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Figure 6. First and second Fourier transforms for 100 squares in 
one field. 



1st Fourier transform 
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2nd Fourier transform 



10 




Figure 5. First and second Fourier transforms for 10 squares in 
one field. 



in Fourier space many rings look different from those pro- 
duced by Gaussian fields. An example of one is shown in 
Figure bl This is t he result of averaging over 10 such fields, 



as in Section 3.3 



this time with a Gaussian background. 
Figure |10| shows the averages for 10 realisations of a Gaus- 
sian random field, for comparison. The characteristic signal 
of the square is still present. 



4.2 Cosmic string maps 

Finally we wish to see if this statistic could detect cosmic 
strings in the CMB. The work is motivated by this, since a 



4 




Figure 7. 100 squares in one field. 



simple approximation for a string map is to consider many 
string segments with random positions and orientations, and 
we hope the ideas illustrated above could apply. We will, 
however, use the results of string simulations to test the 
statistic. 

The simulation gives us fields which have a physical size 
of 2? 75. The Fourier transform of the whole field appears 
Gaussian, so we look at small sections of the simulations to 
see if the statistic picks up the non-Gaussianity. We look at 
fields a quarter of the size of the simulation. We pick them 
at random, and analyse them individually. Figure [H] shows 
two graphs of |-F(/9)| obtained from string maps. Again the 
y-axis is in arbitrary units. Many of the graphs obtained 
look similar to those obtained from Gaussian fields. 

For each string simulation we first obtain \F NG ((3)\, 
then generate 50 Gaussian fields with the same power spec- 
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1st Fourier transform 



Figure 8. 1 square superimposed on a field of Gaussian fluctua- 
tions. 




8x10" 
6x10" 
4x10" 
2x10" 




a 

2nd Fourier transform 
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Figure 9. Average of first and second Fourier transforms for 10 
realisations of 1 square superimposed on a Gaussian background. 



Figure 10. Average of first and second Fourier transforms for 10 
realisations of a Gaussian. 



2nd Fourier transform 



3xl0 6 
2xl0 6 
10 6 




0" 



2nd Fourier transform 




trum as the string simulation and treat each one in ex- 
actly the same way as the string maps to obtain \F^(/3)\, 
i = 1, ...,50. These 50 values of \F? (f3)\ are used to obtain 
the sample variance for each f3: 



where \F° (0}\ = £ £ AFp(0)\. 

Figure n2 shows the number of simulations giving a 

' ' Q 

value of \F((3)\ above or below \F (f3)\ ± cr(f3) for each /?, 
for Gaussian fields and for fields containing strings. The to- 
tal number of string simulations is 20 and the total num- 
ber of Gaussian simulations is 1000. The axes are scaled 
so that the proportion of simulations giving values outside 



Figure 11. |F(/3)| for two different fields containing strings. 



\F (/3)| ± o(0) can be directly compared. Clearly the pro- 
portion is higher for the string simulations, i.e. varies 
more. Large non-zero values of /3 give more variation in 
|-F(/3)| than values close to zero. 

We also want to see if we can detect the strings when 
they are covered in real space by Gaussian fluctuations. We 
have repeated the above procedure for string simulations 
superimposed on a Gaussian field. We have used a power 
spectrum with an exponential fall off on small scales as a 
simple model of the Gaussian fluctuations from before last 
scattering: 1(1 + l)Ci ~ exp(-(Z/1200) 2 ) with the ampli- 
tude of the Gaussian fluctuations being twice that of the 
non-Gaussian fluctuations at I — 500. This power spectrum 
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Figure 12. Number of simulations giving values of |F(/3)| outside 



\F ± cr(f3) for Gaussian fields and fields containing strings. 
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Figure 13. Number of simulations giving values of |-F(/3)| outside 

Q 

\F (f3)\ ± cr(/3) for Gaussian fields and fields containing strings 
superimposed on Gaussian fields. 



can be considered to include noise if it is uncorrelated with 
the signal. This time the results are not so clear. Figure |l3| 
shows the same graphs for string simulations superimposed 
on Gaussian simulations. Comparing the bottom graph to 
the top one would not convince us that it came from a non- 
Gaussian field. 

Analysing the data further, however, we can differenti- 
ate between the two cases. Figure [L4| shows the number of 
simulations giving values of below \F (/3)| — a(f3), for 

Gaussian fields, and for strings and Gaussian superimposed 
fields. For many values of j3, the proportion of |-F(/3)| from 



Figure 14. Number of simulations giving values of \F(J3)\ below 

Q 

\F (fi)\ — <x(/3) for Gaussian fields and fields containing strings 
superimposed on Gaussian fields. 



Q 

non-Gaussian fields below \F (/3)| — a(/3) is higher than the 
proportion that we would expect from Gaussian fields. (This 
is also true for the strings only case, in fact most of the time 
in that case the values of |FX/3)I which are outside the 'one 
sigma' range lie below the lower one rather than above the 
upper one.) 

So to test the sky for non-Gaussianity we need to 
look at several fields, generating Gaussian fields with the 
same power spectrum for each one, and see how many of 
the fields we are interested in give values of |-F(/3)| below 

Q 

\F (j3)\ — <j(f3). Of course we need to look at a reasonable 
number of potentially non-Gaussian fields in order to see if 
the proportion is more than we would expect from Gaussian 
fields. 



5 DISCUSSION 

We have introduced a statistic to look for non-Gaussianity in 
the cosmic microwave background. It is designed to pick out 
non- Gaussian features which are superimposed on Gaussian 
fluctuations and which are therefore not visible in real space. 
Since some scales may be more non-Gaussian than others we 
choose to separate them out. It is o ur experience in previ- 
ous work (Ferreira & Magueijo 1997) that in such situations 
standard statistics, based in real space, fail to recognize the 
non-Gaussianity of the signal. 

Our statistic is naturally tailored for interferometric ex- 
periments, which make measurements in Fourier space. In- 
deed the algorithm exposed above could easily be turned 
into a data analysis package operating over visibilities. 

One of the problems with Fourier space statistics is that 
they recognize only global shapes. Therefore they become 
very ineffective when many individual structures are present. 
By taking another Fourier transform, in angle, over modes 
in a given ring in Fourier space, we can factor out orienta- 
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tions, and be only sensitive to an average shape of individual 
structures. Their random positions, on the other hand, will 
eventually make our statistic very ineffective as the number 
of objects becomes large. We found that although we have 
improved on previous work, still there is a limit on how 
many individual structures there may be in the field before 
the statistic fails to pick out their non-Gaussianity. 

Finally, some comments are in order regarding the ap- 
plications we have considered. The Gaussian component we 
have considered may include noise, if it has the right power 
spectrum. Hence we have already considered the effects of 
a simplified form of noise: it will merely reduce the band 
of scales where the non-Gaussian signal dominates, typi- 
cally providing it with an upper boundary. An exception 
is the case of non uniform noise in the u-v plane, present in 
most interferometer experiments. Non uniform noise in the 
Fourier domain will look non-Gaussian, for it is a Gaussian 
process which is not isotropic or translationally invariant. 
Therefore an extra element of confusion, not dealt with in 
this paper, will appear. 

Regarding the non-Gaussian component we have con- 
sidered two types of signal. We have looked at some of the 
properties of the statistic when it is applied to a field con- 
taining one or more distinct shapes. The geometrical con- 
structions presented are somewhat reminiscent of the jagged 
structures present in small scale dust maps. Experiment- 
ing on fields containing many structures we saw our statis- 
tic rapidly become blind to their non-Gaussianity. However 
interferometer fields are often small. They would therefore 
contain only a few, but more than one, of these structures. 
This is something our statistic can cope with, unlike previ- 
ous Fourier space based statistics. 

We have also studied our statistic when applied to cos- 
mic string maps. We found that some individual small fields 
showed non-Gaussianity in the signal. If the results were 
averaged over many fields the effect of non-Gaussianity on 
\F(/3)\ would not show. Once again this strategy is ideal for 
interferometers, for which combining small fields into a large 
field is in fact a rather awkward operation. 
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